Calculation of parity non-conserving optical rotation in iodine at 1315 nm 
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We examine the feasibility of a parity non-conserving (PNC) optical rotation experiment for the 
2 P3/2 — ^ 2 Pi/2 transition of atomic iodine at 1315 nm. The calculated .EIpnc to Ml amplitude ratio 
is R = 0.80(16) x 10 -8 . We show that very large PNC rotations (greater than 10 /xrad) are obtained 
for iodine-atom column densities of ~ 10 22 cm -2 , which can be produced by increasing the effective 
interaction pathlength by a factor of ~ 10 4 with a high- finesse optical cavity. The simulated signals 
indicate that measurement of the nuclear anapole moment is feasible, and that a 1% PNC precision 
measurement should resolve the inconsistency between previous measurements in Cs and Tl. 
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Introduction - The precise measurement of atomic 
parity nonconservation (PNC) provides a low-energy 
test of the standard model and of internucleon weak 
interactions [TH3] . In recent decades, successful atomic 
PNC measurements have been performed using two 
techniques: (a) the Stark interference technique, on Cs 
U, and (b) the optical rotation technique, on Tl, Bi, and 
Pb [5HH]. The highlight of these efforts was the 0.35% 
precision measurement of nuclear-spin-independent PNC 
in Cs, and the 14% precision measurement of the nuclear 
spin-dependent PNC for the odd-proton nucleus of 133 Cs 
[4]. However, measurement of the anapole moment 
in Cs disagrees with Tl [4j |5], and also with some 
theoretical nuclear calculations ([9HTT] and references 
therein). To help resolve these inconsistencies, and to 
improve the atomic PNC tests of the standard model, 
further experiments are needed. For example, there is a 
proposal to measure the nuclear anapole moment using 
the PNC-generated hyperfine frequency shift of dressed 
states in atoms [12 , with particular proposals for Cs 
[T3] and Fr [14]. However, for other PNC candidates, 
as the precision in the atomic theory is not expected to 
significantly surpass the current experimental or theo- 
retical PNC precision of Cs, future PNC experiments 
have focused on other fruitful directions, such as the 
measurement of atomic PNC on a chain of isotopes 
[9j[15], or the measurement of nuclear anapole moments. 
Therefore, along these lines, PNC experiments are in 
progress on Yb and Dy at Berkeley jTHJ [17], on Rb 
and Fr at the Tri- University Meson Facility (TRIUMF, 
University of British Columbia) [181 HI], on R& + a t 
KVI Groningen [20 , and on Ba + at the University 
of Washington, Seattle [21]. Recently, Bougas et al. 



proposed the measurement of optical rotation in excited 
(metastable) states of Hg and Xe [22], using a novel 
optical bow-tie cavity to enhance available single-pass 
metastable column densities from 10 14 to 10 18 cm -3 
(by achieving ~10 4 cavity passes) and to effect two new 
signal reversals. Hg and Xe both have several stable 
isotopes, and each have two isotopes with odd-neutron 
nuclei (for which anapole moments have not yet been 
measured). It will be useful to find compatible atomic 
PNC candidates that allow measurement of the anapole 
moment of odd-proton nuclei as well, needed to resolve 
the inconsistencies in the Cs and Tl measurements. Here 
we investigate the suitability of ground-state 12T I (which 
has an odd-proton nucleus) as a candidate for measuring 
nuclear spin-dependent PNC. Note that although 127 1 
is the only stable isotope of iodine, several radioactive 
isotopes are also commercially available. Therefore, the 
aim of this paper is to present spin-independent and 
spin-dependent PNC calculations for the 2 P3/2 — ^ 2 Pi/2 
spin-orbit transition of 127 I and to examine the ex- 
perimental feasibility of PNC measurements for this 
transition using the cavity-enhanced optical rotation 
technique. 

Calculations - The Hamiltonian describing parity- 
nonconserving electron-nuclear interaction can be 
written as a sum of the nuclear-spin-independent (SI) 
and the nuclear-spin-dependent (SD) parts (we use 
atomic units: h = \e\ = m e = 1): 

#pnc = #si + #SD 

Gf ( Qw . K j\ ( s /.x 

= T^\~— lb + l aI ) p{rl (1) 
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where Gp ~ 2.2225 x 10 14 a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 

a and 75 are the Dirac matrices, / is the 
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nuclear spin, and p(r) is the nuclear density normalized 
to 1. Within the standard model the weak nuclear charge 
Qw is given by [23] 



Iw 



-0.98777V + 0.0716Z. 



(2) 



Here N is the number of neutrons, and Z is the num- 
ber of protons. The strength of the spin-dependent PNC 
interaction is proportional to the dimensionless constant 
x, which is to be found from the measurements. It is 
believed that k is dominated by the nuclear anapole mo- 
ment [24j [25] . The PNC amplitude of an electric dipole 
transition between states of the same parity \i) and |/) 
is equal to 



(f\d\n)(n\H PNC \i) 



Ei — E n 
(f\H PNC \n)(n\d\i) 



Ef — En 



(3) 



where d = —e J2 i Vi is the electric dipole operator. To 
extract from the measurements the parameter of the 
nuclear spin-dependent weak interaction, x, one needs 
to consider PNC amplitudes between specific hyperfine 
structure components of the initial and final states. The 
detailed expressions for these amplitudes can be found in 
earlier works [26 - 28] . We perform the calculations of the 
PNC amplitude using the configuration interaction (CI) 
method. The PNC amplitude (|3| in iodine is dominated 
by the 5s - 5p transition. Therefore we treat 5s, 5pi/2> 
and 5£>3/2 states as valence states, performing CI calcu- 
lations for seven valence electrons. Calculations are done 
in the V N approximation, which means that the initial 
Hartree-Fock procedure is done for neutral iodine. The 
hole in the bp subshell is taken into account via the frac- 
tional occupation number techniques [29-31 . The com- 
plete set of single-electron orbitals is constructed using 
the B-spline technique [32]. We calculate 40 B splines 
in the cavity of R = 40ae in each partial wave up to 
i max = 2. We use 14 lowest states above the core in 
each of the s, P1/2, <^3/2? an d ^5/2 in the CI cal- 

culations. Core- valence correlations are not included be- 
cause in a many- valence-electron system such as iodine, 
intervalence correlations strongly dominate over the core- 
valence correlations. The effects of external fields, such 
as the electric field of the laser light and electron-nucleus 
parity-nonconserving weak interactions, are taken into 
account within the framework of the random-phase ap- 
proximation (RPA) (see, e.g., [33 ). The summation over 
intermediate many-electron states in (|3| is done with 
the use of the Dalgarno-Lewis technique [34] (see also 
[28} [35]). The main factor contributing to the uncer- 
tainty of present calculations is the incompleteness of the 



many-electron basis. This is very typical for the CI cal- 
culations with a large number of valence electrons. Full 
saturation of the basis would lead to a matrix of unman- 
ageable size. To reduce the size of the matrix, we cut 
our basis in such a way that all important configurations 
are included. First, we limit our single-electron basis to 
l max = 2. Second, we consider only single and double 
excitations when we construct basis configurations from 
the ground-state configuration. Then, we limit single ex- 
citations to 14 lowest states above the core. Usually 14 
out of 40 B splines is enough to saturate single-electron 
excitations. Then we limit double excitations to the eight 
lowest states above the core. Comparing different runs of 
the code with different cuts to the basis we estimate that 
the uncertainty of the present calculations of the PNC 
amplitude is about 20%. This uncertainty can be further 
reduced in case of progress with the measurements by us- 
ing more computer power. The resulting PNC amplitude 
(z component) is 

£1pnc = 1.46(29) x 10- 11 (-Q w /N)iea B . (4) 

Using ([2| and Q we obtain the value of the reduced ma- 
trix element for the £?1pnc transition amplitude between 
states 2 P° /2 and 2 P° /2 of 127 I: 



£1pnc = 3.35(67) x 10- n iea B . 



(5) 



The spin-dependent PNC amplitudes between different 
hyperfine-structure (hfs) components of the ground and 
excited states are presented in Table [T] The Ml ampli- 
tude between the ^3/2 an d 2 P?/ 2 states is 0.00420 ea B . 
Here we took into account that the overlap of the 
5pi/2 and 5£>3/2 functions is found from the Hartree- 
Fock calculations to be (5pi/ 2 | 5£>3/2) =0.997. The 
uncertainty for the Ml amplitude is order 0.1%. The 
angle of the optical rotation is proportional to the ratio 
R = Im(£pNc/Ml). Using ^ and Ml = 0.0042ea B 
we get R = 0.80(16) x 10" 8 . Note that the dependence 
of the £pnc and Ml amplitudes on the values of the 
total angular momentum F (including nuclear spin) 
and its projection M is the same if we neglect the k 
contribution to £Wc and the nuclear magnetic moment 
contribution to Ml. Therefore, for the nuclear-spin- 
independent contributions the values of R are the same 
for all transitions in Table [U 



Experimental feasibility 
is given by 



PNC optical rotation 



^pnc = -~y [ n M -i\R 



(6) 



where I is the length of vapor, A is the optical wavelength, 
uj is the optical frequency, and n(uj) is the refractive in- 
dex due to the absorption line. For optical depths L <C 
1, the optical rotation angle (^pnc pi (where p is the 
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TABLE I: PNC amplitudes (z components) for the 



|5p 5 



3/25 



Fl , Mrr 



minimi, i 7 ^}) transitions in 127 I. .F — J+J, where J is nuclear 
spin, / = 5/2. 



PNC amplitude [10" n zea B ] 



0.92(18) [1 -0.043(9)x] 
1.21(24) [1 -0.045(9)x] 
0.39(8) [1 + 0.037(7) x] 

-0.64(13)[l-0.048(10)x] 
1.08(22)[1 + 0.035(7)x] 

-0.84(17)[l + 0.031(6)x] 



density and pi is column density of the vapor). For op- 
tical depths L ^> 1, the vapor is optically thick near the 
line center where ^pnc is largest and can no longer be 
observed. The effective maximal rotation angle is shifted 
further off resonance as \fpl, and ^p^c ^ V7^? as can 
be shown by maximizing the product of dispersion and 
transmission. Therefore, the rotation angle can still be 
increased with increasing column density for L ^> 1. For 
this reason, in the PNC optical rotation experiments us- 
ing Tl, Bi, and Pb vapors, optical depths of about 50 
were used (corresponding to column densities of ~ 10 19 
cm -2 ) to obtain optical rotations of ~1 /xrad. The value 
of R = 0.8 x 10 -8 for 127 I is about an order of magnitude 
smaller than that for Tl, Bi, and Pb (whereas the Ml am- 
plitude and transition wavelength A are about the same in 
all four cases, and the line-shape parameters are approxi- 
mately similar). Therefore to compensate for the smaller 
value of R and to produce PNC optical rotations of ~1 
/xrad for 127 I, larger optical depths of ~ 10 3 are necessary 
(corresponding to column densities of order ~ 3 x 10 20 
cm -2 ). Such a large column density can realistically be 
achieved only with the use a high-finesse cavity. We first 
discuss the feasibility of producing column densities up 
to ~10 22 cm -2 , and then use such column densities in 
simulations of optical rotation profiles for 127 I. 



High iodine atom densities of ~ 10 cm -3 have been 
achieved in dc glow discharges [36-38]. However, in 
achieving these densities through discharge, relatively 
high pressures (~few tens of Torr) of both the precur- 
sor (I2, CH3I) and carrier gases (He, Ar) are required. 
Such pressures can significantly suppress the measured 
PNC signal, through the induced collisional broadening 
of the 2 P3/2 — ^ 2 Pi/2 transition. Therefore, we consider 
an alternative method for the creation of ground-state 
2 P 3 / 2 iodine atoms at lower pressures through photodis- 
sociation of I2 molecules: 



I 2 



I( 2 P3/2)+I( 2 P 3 /2) 



The absorption cross section of I2 at 308 K peaks at 
a = 2.8 x 10" 18 cm 2 for A -530 nm (see Fig. [l]and [39]). 
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FIG. 1: (A) Photodissociation of I2 (ideally with light at 
^532 nm) and the corresponding absorption cross section 



(taken from [39]). (B) Level scheme for the P 



3/2 



/2 



Ml transition of atomic iodine at 1315 nm. Indicated in gray 
color are the individual F — »• F' transitions constituting the 
two separated hyperfine groups of Fig. [2] The main contri- 
bution to PNC mixing comes from the 5S core states. 



The production of spin-orbit-excited I( 2 P!/ 2 ) atoms: 

I 2 ^>I( 2 P 3/2 )+I( 2 P 1/2 ), (8) 

is energetically allowed only for photodissociation wave- 
lengths smaller than ~500 nm, so this undesirable chan- 
nel is avoided by photodissociating the I2 molecules at 
higher wavelengths, e.g., at 532 nm, with the second har- 
monic of a Nd:YAG laser or with a green light-emitting 
diode array. At room temperature, the vapor pressure of 
solid I2 yields a density of about 5 x 10 15 cm -3 I2 vapor 
and a column density (over a 100-cm path length) of 5 
x 10 1T cm -2 . This corresponds to an optical depth of 
~1. The I( 2 P 3 / 2 ) atoms recombine with the three-body 
reaction: 



I( 2 P 3 / 2 )+I( 2 P3/2)+l2 



2h 



(9) 



where the recombination rate is given by k r = 3.1(3) x 
lO -30 cm 2 molecule- 2 s -1 at 308 K [40 . Setting the 
photodissociation rate and recombination rate of I2 to 
be equal, we obtain: 



d[h 



dt 



= -J[I 2 ]+fc r [I 2 ][I] 2 =0, 



(10) 



where J is the I2 photodissociation rate constant, given 
by J = $(7, and is the photon flux. Solving for the 
atomic iodine concentration yields [I] = yj &a/k r . Choos- 

1 (corresponding to 50 W 



ing <E> ~ 10 2U photons cm s" 

cm -2 at 532 nm, which is probably close to a reasonable 
(7) upper bound) gives a steady-state iodine- atom concen- 



tration [I] ~ 10 16 cm -3 . At I densities of 10 16 cm -3 , the 
mean time between wall collisions (for a cell-tube diame- 
ter of about 1 cm) is ~1 ms. The sticking probability of 
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I atoms with the cell walls is typically on the order of 5% 
but can be reduced to 10 -3 for acid-coated cells [4T| |42], 
so that loss of I atoms from sticking and recombination 
at the walls should occur on the ~l-s time scale. The I 
atoms should collide with the cell walls hundreds of times 
on average, ensuring that they are thermalized. Notice 
that we calculate a steady-state concentration [I 2 ] =0 
(the initial [l2] = 5 x 10 15 cm -3 has been photodissoci- 
ated to give [I] = 10 16 cm -3 ), as we have ignored recom- 
bination processes with smaller rates in Eq. (10), such 



as I-atom recombination at the walls and three-I-atom 
recombination. However, as these reactions are expected 
to be less significant than ([9| at these pressures, the fi- 
nal I2 concentration can be kept at ~10 15 cm -3 . Finally, 
using a cell with a length of 100 cm and a high-finesse 
optical cavity with ~10 4 passes gives an upper bound for 
the effective I-atom vapor column density of 10 22 cm -2 . 

In a typical PNC optical rotation experiment, the ob- 
servable is the product of (^pnc times the transmitted op- 
tical power and this is doubled by the signal reversals pro- 
posed in [22 . Simulations of PNC optical rotation in the 
2 P3/2 ^ 2 Pi/2 ^1 transition in iodine are shown in Fig. 
[2] In Fig. [2|a) we assume a column density of 10 22 cm -2 
and the dependence on the Lorentzian broadening of the 
line is investigated. We see that an increasing Lorentz 
width reduces the measured PNC signal. As the Lorentz 
width is mainly determined by the pressure inside the 
cell, maximizing the ratio of iodine-atom pressure to to- 
tal pressure is desirable. The collision-broadened Lorentz 
widths for the Ml transition for I and I2 collision part- 
ners were estimated in [43: to be 4.3 and 10.3 MHz/Torr 
at room temperature, respectively. We see that, for re- 
alistic Lorentz widths, signals in excess of 10 /irad can 
be achieved, along with sufficient hyperfine resolution, 
which is necessary in determining the x- value of Table [l| 

In Fig. |2jb) the reverse situation is examined: a con- 
stant Lorentz width is maintained, determined by fixed 
densities of I and I2, and a varying number of passes is 
assumed, ranging from 10 4 down to 300. As we can see, 
even in the case of 300 passes, a signal on the order of ~1 
/irad at the peak is predicted (similar to the signal levels 
measured for Tl [5]), allowing for a significant margin of 
error in the experimental constraints. 

Finally, in Fig. [2|c), the nuclear spin- dependent PNC 
effects are examined. We use the previously measured 
values of k from the thallium [5] and cesium [4] experi- 
ments to predict values for 127 1. Using the simple single- 
valence-proton model [25], we obtain: 



-( 127 I) 
x( 127 I) 



x( 205 Tl) 
- x( 133 Cs) 



0.13(20) [from Tl] (11) 
- 0.38(6) [from Cs] (12) 



These two values define the anticipated physical range of 
the spin-dependent part of the measurement and are the 
ones used for the simulations. We see that peak signal 
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FIG. 2: Simulated signal vs frequency, (a) For the same 
column density (10 22 cm -2 ) and various Lorentz contribu- 
tions (from 1 to 50 MHz), (b) For the same Lorentzian width 
(3 MHz) and various column densities (from 10 22 down to 
3xl0 20 cm -2 ), (c) Effect of the nuclear anapole moment. 
This is the simulated signal for p = 3 x 10 21 cm -2 and Tl =3 
MHz, for k — 0.13 and h — —0.38. For a discussion, see text. 



values differ by about ±2-3% for the two cases, result- 
ing in total signal differences of up to ~5%, as expected 
from the h coefficients in Table |TJ Therefore an experi- 
mental precision of 1% or better is needed to clearly dis- 
tinguish between these two extremes. Also evident from 
Table [i] and reflected in Fig. |2|c) is the fact that tran- 
sitions belonging to the two separated hyperfine groups 
[F -» F' = 2 and F -» F' = 3) deviate in opposite 
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directions in the two cases, meaning that a nuclear spin- 
dependent signal could be acquired even in the case where 
broadening mechanisms would blur the ground-state hy- 
perfine structure. 

We note that within the same single-valence-proton 
model, x( 129 I) ~ x( 131 I) ~ x( 133 Cs), so that measure- 
ments on these iodine isotopes could provide additional 
experimental tests. The nuclear magnetic moments of 
these three nuclei (all with nuclear spin I = 1 12) agree 
within 2% [44-46 , which supports the validity of the 
single- valence-proton model. 

Conclusions - Calculations for atomic iodine 
presented here along with signal simulations, indicate 
that a PNC optical rotation experiment is feasible, 
utilizing the cavity-enhanced, optical-rotation scheme 
of [22]. A nuclear anapole moment measurement seems 
possible, as the simulated signal levels and hyperfine 
resolution (shown in Fig. [2| and the calculated ampli- 
tudes of nuclear spin-dependent PNC (shown in Table 
[i]) compare favorably to those of the successful PNC 
optical-rotation experiment in thallium [5] . 
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